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Abstract 



We prove the existence of a large class of dynamical solutions to the Einstein-Euler equations that have 
a first post-Newtonian expansion. The results here are based on the elliptic-hyperbolic formulation 
of the Einstein-Euler equations used in [15], which contains a singular parameter e — vt/c, where 
■ ut is a characteristic velocity associated with the fluid and c is the speed of light. As in [15], energy 

estimates on weighted Sobolev spaces are used to analyze the behavior of solutions to the Einstein- 
Euler equations in the limit e \ 0, and to demonstrate the validity of the first post-Newtonian 
expansion as an approximation. 



1 Introduction 

The Einstcin-Euler equations, which govern a gravitating perfect fluid, are given by 

> , 

5^ ■ G'^ = ^T*^ and V.r^ = 0, 

l> 

CO ' where 

T'^ = (p + c-^pyv^ +pg'K 

with p the fluid density, p the fluid pressure, v the fluid four-velocity normaUzed by v^Vi = — c^, c the 



, speed of Ught, and G the Newtonian gravitational constant. Defining 

Vt 

, e = — 

X, 

^ , where vt is a typical speed associated with the fluid, the Einstein-Euler equations, upon suitable rescaling 

[15], can be written in the form 

G-y = 2e^T'^ and V,r^ = 0, (1.1) 

where 

T'^ = {p + e^pWv^ + pg"-^ and w'wi = 

In this formulation, the fluid four- velocity v^, the fluid density p, the fluid pressure p, the metric gij, 
and the coordinates (x*) i = 1,...,4 are dimensionless. By assumption, the (a;*) are global Cartesian 
coordinates on spacetime Af = M'^ x [0, T), where the (x^) (/ = 1, 2, 3) are spatial coordinates that cover 
M."^ , and t = x^/vt is a Newtonian time coordinate that covers the interval [0, T). By a choice of units, 
we can and will set ut = 1- 
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Post-Newtonian expansions for the Einstein-Euler system refer to expansions of solutions to this 
system in the parameter e, about e = 0, where the lowest expansion term is governed by the Poisson- 
Euler equations of Newtonian gravity: 

dtl + di{p^') = Q, (1.2) 

"pidtw^ + w^diw'') = -{pd-^l? + d^p) , (1.3) 











A$ = p. (1.4) 

Here p, p, and w'^ are the fluid density, pressure, and three velocity, respectively. 

Formal calculational schemes for determining the post-Newtonian expansion coefficients and the equa- 
tions they satisfy exist, and are in wide use by physicists [5,9]. In fact, these post-Newtonian compu- 
tational schemes are one of the most important techniques in general relativity for calculating physical 
quantities for the purpose of comparing theory with experiment. For example, in gravitational wave 
astronomy, post-Newtonian expansions are used to calculate gravitational wave forms that arc emitted 
during gravitational collapse [5] . 

It is important to stress that the formal post-Newtonian expansion schemes all implicitly rely on 
the assumption that the expansions exist and approximate solutions to general relativity. Therefore, 
to establish existence of such approximations, and to answer questions about their range of validity, a 
different approach must be taken to the problem. In [15], we took a first step in analyzing this problem 
by proving the existence of a wide class of one-parameter families of solutions to the Einstcin-Euler 
equations that converged in a suitable sense to the Poisson-Euler equations in the limit e \ 0. We also 
remark that similar results were also established, using a different method, by Alan Rendall [19] for the 
Einstein- Vlasov equations. 

In this paper, we use the results of [15] to prove the existence of a large class of solutions to the 
Einstein-Eulcr equations that can be expanded in e to the first post-Newtonian order. Moreover, we 
demonstrate the existence of convergent expansions in e for solutions to the Einstein-Euler equations. 
These expansions are, in general, not of the post-Newtonian type since the expansion coefficients can 
depend on e. Nevertheless, the expansions are convergent, and therefore, represent a kind of generalized 
post-Newtonian expansion. We note that analogous expansions for the Vlasov-Maxwell equations and 
Vlasov- Nordstom equations have been rigorously analyzed in [2-4] . 

The difficulty in analyzing the post-Newtonian expansions arise from the fact that the limit e \ is 
singular. To analyze this limit, we follow the approach of [15], which requires that the metric Qij and 
the fluid velocity are replaced with new variables that are compatible with the limit e \ 0. The new 
gravitational variable is a density u*-' defined via the formula 

where 

From this, it not difficult to see that the density u^^ is equivalent to the metric Qij for e > 0, and is well 
defined at e = 0. For the fluid, a new velocity variable is defined by 



^ and w^ = - -. (1.7) 



V = w 

e 



For technical reasons, we assume an isentropic equation of state 

p = JTp^^+i)/", (1. 
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where K e IR.>o, n € N. This allows us to use a technique of Makino [14] to regularize the fluid equations 
by the use of the fluid density variable 



{AKn{n+l)) 



a^". (1.9) 
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The resulting system can be put into a symmetric hyperbolic system that is regular across the fluid- 
vacuum interface. In this way, it is possible to construct solutions to the Einstein-Euler equations that 
represent compact gravitating fluid bodies (i.e. stars) both in the Newtonian and relativistic setting 
[14,18]. In the Newtonian setting, this is straightforward to see. Using (|1.8|) and p.9p . the Poisson-Euler 
equations (|1.2p - (|1.4p imply that 



(1.10) 

(1.11) 
(1.12) 

which is readily seen to be regular even across regions where a vanishes. 

As discussed by Rendall [18], the type of fluid solutions obtained by the Makino method have freely 
falling boundaries and hence do not include static stars of finite radius, and consequently this method 
is far from ideal. However, in trying to understand the post-Newtonian expansions, these solutions are 
general enough to obtain a comprehensive understanding of the mathematical issues involved in the 
post-Newtonian expansions. 

As in [15], our approach to the problem of post-Newtonian expansions is to use the gravitational 
and matter variables {u*^ , w', a} along with a harmonic gauge to put the Einstein-Euler equations into a 
singular (non-local) symmetric hyperbolic system of the form 

\P{tW)dtW = \:^diW + h\e, W)diW -f i^(e, W). (1.13) 

Singular hyperbolic systems of this form have been extensively studied in the articles [6, 11, 12, 20, 21]. 
Especially relevant for our purposes, is the paper [21]. There, a systematic procedure for constructing 
rigorous expansions to singular symmetric hyperbolic systems is developed (see also [11,12]). However, 
the techniques of [6, 11, 12,20,21] cannot be applied directly to our case. The reason for this is that the 
initial data for the system (jl.l3p must include a 1/r piece for the metric and cannot lie in the Sobolev 
space li^ . This problem was overcome in [15] by using a one parameter family of weighted Sobolev 
spaces that include 1/r type fall off for e > 0, and reduce to the standard Sobolev spaces in the limit 
e \ 0. We again use these weighted Sobolev spaces, this time to generalize the results of [21] so that we 
can apply them to the problem of generating rigorous post-Newtonian expansions. 

The next theorem is the main result of this paper, and the proof can be found in section [51 The 
definition of the spaces and XT,s.k.5 can be found in Appendices \K\ and [B] 

Theorem 1.1. Suppose -I < 5 < -1/2, s > 3, fc > 3 + s, a,w\i{'' e H^_,, f e H^Io ' suppa C Br, 

o o o 

d let Tq ^ is the maximal existence time (see Proposition \3. ?| ) for solutions to the Poisson-Euler-Makino 



ana 



equations (|1.10p - p.l2p with initial data a(0) = a, {{)) = . Then for any Tq < T there exists an 
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eo > 0, and maps 

ul^{t) : ul^{t)^u'^{0), diul^it), dtul'{t)&XT„,s.k.5-i 0<e<eo, 
a^{t), wl{t) e XTo,s^kJ-i < e < eo, 

r 

a{t), W {t) e XTo.s,k,S-l, ^{t) £ XTo,s,k+2J with dt^{t) e XTo,s,k+l,S-l, 
I'^it) : b{t) - I'^iO), dll'^it) e XTo.s-g.k-g.S-l 9 = 1, 2, 
a{t), W\t) £ XTo,s-q.k-q.S-l <7 = 1,2, 

k^{t) : k'it) - u:^ (0), dik^it) e Xto^s-3M-3,5-i {q, e) E Z>3 x (0, eg], 
a^{t), wl{t) e XTo,s-3,k-3,s-i {q,e) e Z>3 x (0,eo], 

such that 

(i) the triple {u*-' (i), Q;e(i), determines, via formulas (ll.5p - (|1.9p . a solution to the Einstein- Euler 
equations p.ip m f/ie harmonic gauge for < e < eo on the spacetime region (x^ ,t = x'^) E D = 

R3 X [o,ro), 

(ii) dtu'/{0) = e^ii-\ dfv.i-^{0) = e^f", a,{0) = a, and (0) = for < e < eo, 

o o 

" 

(Hi) {a{t),w'{t),<^{t)} is the unique solution to the Poisson-Euler-Makino equations (|1.10p - p.l2p with 
initial data a(0) = oi, (0) — , 

o o 

(iv) forq = 1,2, {\x'^ {t),a{t),w\t)} satisfi es a linear (non-local) symmetric hyperbolic system that only 

00 ^ 00 ^"^- ll- 

depends on {a{t),w\t),^{t)} ifq = l, and {a{t),w^ (t), ^{t),VL^^ {t),a{t),w''{t)} ifq = 2, 
(v) for q G Z>3, {u^J (t) , ag(t) , wl(t)} satisfies a linear (non-local) symmetric hyperbolic system that 

OOr^ P ■ P P ■ P ■ ■ P P ■ 

only depends on e, {a{t),w {t),^{t)}, {u^^ (t) , a{t) , {t)} for p ~ 1,2, and {u^^^ {t), ai{t),wl{t)} 
/or p = 3,4, ... ,g - 1, 

(vi) {u^^ , aii{t),wl{t)} and {iil^ (t), ae(t),wl{t)} for q G Z>3, satisfy the following estimates: 

\\uHt)\\L. + ||9,u:^(t)||^. +e||9iu:-'(i)||^. +e||aAu:^^^ < 1, 

\\a,{t)\\H>^ + \\wl{t)\\H^ + \\dta,{t)\\H>^-i + \\dtwl{t)\\H>^-i < 1, 

wkmLi + ||a/u:^'WIU-3 + e||ai:^"(0llif.-3 +e||5,azo:^(t)ll^.^^ < 1, 

m-^ ^ 1: 

for all (i, e) G [0,To) x (0,eo], and 
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(vii) {u^^ , a^{t) , wl{t)} admits convergent expansions (uniform for < e < eo) of the form 



9=1 9=3 



e 1 



P 2 oo 

q=l q=3 
^ oo 

q=l q=3 
2 oo 

q=l q=3 
2 oo 

where the first expansion is convergent in C'^{[0,To); Lg), and the rest are convergent in both 
CO([0,ro);i/'^-4) andC"{[0,Toy,H^:lJ. 

Remark 1.2. 

(a) For q — 1,2, the equations satisfied by {u^^ ,a,w'^} are the ones obtained by directly substituting 



the expansions of Theorem 1 1.11 (vii) into the Einstein-Euler equations and collecting terms to order 
e^, and therefore coincide with the standard first post-Newtonian expansions. 

(b) The equations satisfied by {iXe^^ ,ae,We^} for q> 3 can be determined from the equations satisfied 

9 

by the We defined in the proof of Theorem 15.11 

To facilitate comparisons of the approach taken in this paper with previous studies, we define the 
following e-independent quantities: 

where (rjij) = diag (1, 1, 1, —1). Then a straightforward calculation, using statement (vii) of Theorem ll.il 
and formulas (|1.5p - (|1.6p . shows that the metric gtj can be expanded as follows 



544 - - 2$ - eb^ - (3(1)' + b^) + O(e^), 



and 



gij = Sij - 2e'Sijl> - e'h" ~ e^[{l>fSij + P'') + 0(6^). 

It is worthwhile to note that higher order expansions in e can be generated for the metric gij using part 
(vii) of Theorem ll.il These higher order terms will, in general, depend on e in a non-analytic fashion, and 
therefore, without further analysis, the relation of these expansion terms to the standard post-Newtonian 
expansions is not clear. 



2 Einstein-Euler equations 

In this section, we quickly review the formulation of the Einstein-Euler equation used in [15] to analyze 
the limit as e \ 0. 
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2.1 Reduced Einstein equations 

As discussed in the introduction, we use a symmetric tensor density u*-' instead of the metric , which 
for e > completely determines the metric via the formula 



1 



where 
and 

To fix the gauge, we let 
and demand that 



gij ^ij + 4g2yy ^ _ _ det(r^). 



dk 



I3X3 

-1 



di if fc = / 
edt if A: = 4 



diU'^ = 0. 

For e > 0, this condition is easily seen to be equivalent to the harmonic gauge 

dkd"^ = 0. 



Here g^^ = \J— det{gke)g^-' is the metric density in the coordinates (a;'). 
Next, defining 



u 



^3 



(flu) 



and 



(2.1) 
(2.2) 



(2.3) 

(2.4) 



(2.5) 
(2.6) 
(2.7) 
(2.8) 
(2.9) 
(2.10) 



(2.11) 
(2.12) 



the Einstein equations G'^ = 2e*T'-' , in the harmonic gauge, can be written in first order form as 



(2.13) 
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where 



and 



with 



/I - 4eu44 0\ 

A\eu)^[ 6'-' + Aeu'-' , (2.14) 

V 1/ 

/ S^'^ 0\ 

= U^'^ , (2.15) 

V 0/ 

/4u« 4u^-' 0\ 

A^(u) = 4u" , (2.16) 
\ 0/ 

f;^(u) = (0,0,<^")^, (2.17) 

FY (u, eu) = (A*^' + B*^' + C'^ , 0, 0)^, (2.18) 



JiJ+s'' (2.19) 



e 



e-'p 



^^^'^ = P [\g\vJv* [(|.9| - l){v^r + {{v^f 1)] 

+ 1^ tpv'v^ + 4e|.gri/2pu4j ^(^4)2 ^ |g|-i/2p(_i + 4£u44) J ■ ^^-^^l 

Letting 

w=(a,w')^, (2.21) 

we can decompose S^^ as 

S''^Sl'+eS\\ (2.22) 

where 

5^^(u,w,eu,ew)=p(^l_l^,,^^°^^^4^ e"! [(|g| - l)(l' W^)2^+7(l + 6z«Y " l)]) ' ^^'^^^ 

and 

<?y7w .wA - /^^'"^'"■^+^^"'^^^^ + pew^(l + ew'')+4|5|-i/2peu^4 \ 

^,1 (w, eu, ewj - ^^^j ^ ^^^^ ^ A\g\-^/V^ p(l + ew^ + \g\-^^M-^ + 4eu44) J ' ^^■^^> 

We will refer to the gauge fixed Einstein equation (j2.13p as the reduced Einstein equations. Because of 
the matrix inversion (|2.9p used to define the inverse density Qij, the reduced Einstein equations will be 
well defined provided 

eu G V = { (r'J ) e M4X4 | det(77*^' + Ar'^) > } . 

2.2 Euler equations 

In [15], we also showed that if we use the fluid variables (|2.2ip . and choose initial data that satisfies 

= := evy + 1/e = eg44(l/e + w")^ + 1/e + 2g4j(l + ew^)w' + egijw^w\ (2.25) 
then the Euler equations V iT^^ = are equivalent to the system 

a*94W = a^div^ + b, (2.26) 
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where 



r^ = ^r, i9^,)^i9''r\ (2.28) 

L] = (5; + e^v'vj , = (2.30) 

Mij=g,j+2e^v,vj, (2.31) 

f y - (fl'"(20.^fe - 0z,0£p)5™u^P + 2{Qip5\^d^)u'P - 20,(,9,)u'=^)) , (2.32) 



and 



" =( eqLf M,,{l + ew'))^ ^^.33) 



b-{ 'fM^A- (2.35) 



We also note that 

>0 (5, 



a^=[\ j +a^(eu,ew), (2.36) 



a 



a_ci _c ■'i ] + a{eu, ew) + aa^ {eu, ew) , (2.37) 



and 



\ + ( '^^i('^"'^'*^) ■ '^"fc 

-77"" {2i-iu-qip + ?7£p) uJP - 2 (?7£p(5|uf - 2r)i4uf) J \ b2 (eu, ew) . Ufc 



(2.38) 



where {a^, a, , bi, 62} are analytic in all their variables provided that eu G V, {a^, a, a^} are symmetric, 
and d^{0, 0) 0, a^(0, 0) = 0, a(0, 0) = 0, 61 (0, 0) = 0, and 63(0, 0) = 0. 

3 Uniform existence and the zeroth order equations 

The combined systems (|2.13p and (|2.26p can be written as 
b°{€V, €^U)dtV = -cJdiV+b'iV, eU, eV, e^U)diV + foiV, eU, eV, e^U) + efi{V, eU, eV, e^U) + -giV), (3.1) 
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where 



C/=(0,0,u^^O,0)^ u'i=u'i\ (3.2) 

o o ' 

V = (u4^uj', (Su*^ a, w'f , Su'^ = u'^ - eu*^ (3.3) 



and 



a (eu, ew) 

^^-fo' o)' (3-^) 

foiV^eU,eV,e^U) ^ pQ-' W " w, eu, 6w)\ 

■' ^ ^ \ 6(u,w, eu, ew) / ' ^ ^ 

/i(y,ef/,ey,e2C/) = Min, en) - Sj , eu, ev.)\ ^3 



5(^) = (-<5^<5ip(a),0,...,0)^. (3.9) 



For initial data, we will often use the following notation: given a function z that depends on time t, we 
define 

z = z\t=Q . 

o 

In addition to solving these equations, we must also solve constraint equations on the initial data to 
get a full solution to the Einstein-Euler equations. Letting 

= g'^^Bleu'^ + e^ {A'^ + B'^ + C'^) + q'^BI^u*''^ - 2a^^u'=(*0^)^ (3.10) 

and defining 

C'^ = e-i(g4'^-T4'0, c4 = g44_^44^ W^d.u'', 
the constraint equations to be solved on the initial hypcrsurfacc Sq = {{x^ , 0) | (x^) S R^} are: 

= (gravitational constraint equations), (3-11) 
W = (harmonic gauge condition), (3-12) 



and 



TV = (fluid velocity normalization). (3.13) 



To fix a region on which the system where both the evolution p.ip and constraint equations p.lip - 
p.l3p are well defined, we note from (|2.14p . (|2.36p . and the invertibility of the Lorentz metric {rj''^) that 
there exists a constant Kq > such that 

-det(77*^' +4eu^^) > 1/16, l + ew^>l/16, (3.14) 

A\eu)>^l, a4(eu,ew)>-^I, (3.15) 
16 16 

and 

|A4(eu)|<16, |a''(eu,ew)| < 16 (3.16) 

for all |eu| < 2Ko, |ew*| < 2Ko, \ea\ < 2Ko. The choice of the bounds 1/16 and 16 is somewhat arbitrary, 
and they can be replaced by any number of the form 1 /M and M for any M > 1 without changing any 
of the arguments presented in the following sections. However, since we are interested in the limit e \ 0, 
we lose nothing by assuming M = 16. 
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3.1 Newtonian initial data 

In [15], we proved the following theorem, based on previous work by Lottermoser [13], concerning the 
existence of e-analytic solutions to the constraints (|3.1ip - (|3.13p . Before we state the theorem, we note 
from (|1.9p . (jl.Sp . and the weighted multiplication inequality (see [15] Lemma A. 8 ) that if a e Hg 
i5<0,k> 3/2) then p,pe H^. 

Proposition 3.1. Suppose -1 < (5 < 0, fc > 3/2 + 1, i? > and {p,p,w^ ,j" € {H^Z^? x H^-i x 
HgZi X Bfi{Hf) . Then there exists an eo > 0, an open neighborhood U of {p,p^w^ ,1^'^), and analytic 
maps {-eo,eo) X U Hl_^ : {t, p,p,w' .i^;' w'^ , {-eo.eo) x U -> : {e, p,p,w' 
(f>, (-eo,eo) X U -^H^ : {e, p,p,w' such that for each {p,p,w' ,}^'') £ U, 
{e, p,p,w^ jW^jU^ jd^u^^) is a solution to the three constraints 



whe 



^0, W = 0, and TV 0, (3.17) 



e £544 
Moreover, if we let (po — 4>\e=o, Ifo = rD^|e=Oj andwQ — w*\i=o, then (po, ItIq, andwQ satisfy the equations 

A(/)o = p, Arog = -(9l34'^ + pw^, ^^0 = 0, 

respectively. 

In section [SJ we show that the analytic dependence of the initial data on e implies that there exists a 
corresponding convergent expansion in e for the solution generated from the initial data. 

3.2 Uniform existence 

To prove local existence of solutions to ()3.ip on a uniform time interval independent of e, we use a 
non-local symmetric hyperbolic version of p.l|) . This system is essentially the one used in [15] to derive 
uniform existence, convergence, and error estimates for the limit e \ of solutions to p.ip . However, 
we employ a few refinements that can be used to simplify the proof in [15], and will also be useful for 
analyzing the higher order expansions in e. 

Letting Xr G Cq" be a cutoff function that satisfies 

Xr\b^= 1, < Xii < 1 , and suppXi? C B^r , 

we replace g{V) in (|3.ip with 

g{V) = i-5l6ixRpia),0,...,0f , (3.21) 
and, following [15], we define the Newtonian potential by 

A$ = XflP. (3.22) 
Before proceeding, we first recall the following inequalities from [15]: 
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(a) If ^ > 3/2, there exists a constant Csob such that 

ll-lli- <CsobM|H£, Vee[0,eo]. (3.23) 



(b) For eo > and < -3/2, 

(c) For eo > 0, and -2 < r] < -3/2, 

(d) For eo > and > -3/2, 
(d) If £2 < £1, and 771 < ?72, then 



k,^,., < II ■ Ik^ VeG[0,eo]. (3.24) 

U^<I|-|U^„ Vee[0,eo]. (3.25) 
Vee[0,eo]. (3.26) 



(3.27) 



Lemma 3.2. Suppose eo > 0, — 1 < 77 < —1/2, and £ > 3/2. Then the maps 
and 

a, o $ : — . : a ^ a,$(a) 

are uniformly analyti^ for e G [0, eo]. 

Proof. First we recall that for —1 < 77 < —1/2, the Laplacian 

A : i7f 2 ^ (3.28) 

is an isomorphism by Proposition 2.2 of [1]. Next, by assumption £ > 3/2, and hence it follows that the 
map J 9 a ^-> p = {4:Kn{n + l))^"a^" e ^i^-i c is uniformly analytic for e G [0, eo] by Lemma [A. 71 

Moreover, the linear map H^_i ^ 9 it 1— > Xr''^ ^ ^^-2 is clearly well defined and uniformly bounded for 
e e [0, eg]. Since compositions of uniformly analytic maps are again uniformly analytic, we see that the 
map 3 a A^^(x^/3(a)) G H^'^'^ is uniformly analytic of e € [0,eo]. 

Next, we recall that differentiation H^~^^ 9 u 1— > dju e H^^'^\ is a bounded linear map, and the 
imbedding H^t.i C H^'^\ ^ is well defined and uniformly bounded for e g [0,eo] by p.24p . Again using 
the fact that uniform analyticity is preserved under compositions, we get that the map 9/ o <i> : H^_i^^ 
^ri-i.e is uniformly analytic for e G [0, eo]. □ 

Following [15], we use the Newtonian potential to define a new combined gravitational-matter variable 
W via the formula 

W^V-d'^?, (3.29) 

where 

d$ := (0, SlSidj'i>{a), 0, 0, 0) . (3.30) 

Notice that the transformation (|3.29p leaves the matter variables unaffected. Consequently, we can define 
W by 

and treat $ or d^ as a function of VF. In fact, by Lemma [3T2l 



hLi., 3W^d<i>e (3.31) 



^See Appendix lAl for a definition of the term uniformly analytic 
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defines a uniformly analytic map for e e [0, eo]. 

To formulate the evolution equation entirely in terms of W, we need the "time derivative" of the $ 
map. So we define 

<^{W, eU, eW, e^U) A'^ f __^g___n(a4(eu, ew)-i [a^(w, eu, ew)diw + 6(u, w, eu, ew)] ) 

(3.32) 

where n((a,w')"^) = a is a constant projection map. By construction, $ = dt^ when evaluated on a 
solution of the reduced Einstein-Euler equations. 

Lemma 3.3. Suppose Ri > 0, eq > 0, —1 < i] < 1/2, and £ > 3/2. Then there exists an R2 > such 
that the maps 

$ : Bn.iH^^^^J X BR.iH'^) X Br,{H'^^_,^^) X BR.iH^) H'^+' : {W,U,W,U) ^ Mw,U,W,U) 
and 

5/0$ : Br,{H'^_,JxBr,{H'^)xBr,{H'^_,^^)xBr,{H'^) H'^^,^^ : {W,U,W,U) ^ di{MW,U,W,U)) 
are uniformly analytic for e £ [0, eq] • 

Proof. Fixing i?i > 0, eg > 0, -1 < 77 < -1/2 and £ > 3/2, it follows directly from Lemmas [X2l and [Ol 
that there exists a i?2 > such that the map 

Br,{H^,_,^,) X Br,{H'^) X Br,{H'^_,J X Br,{H^,) 3 {W,U,eW,e^U) ^ 

X7ja^""^n(a''(eu, ew)"-^ [a^(w, eu, ew)diw + 6(u, w, eu, ev^f)] e H^zl 

is uniformly analytic for e G [0, eg]. The rest of the proof now follows from the same arguments used in 
the proof of Lemma [3?2l □ 



To fit with the above notation, we define 

d$ = (0,4(5^(9/$, 0,0,0)^. 

Noting that 

b°{eV,e^U) ^b°{eW,e^U) and b' (V, eU, eV, e^) ^ {W,eU,eW, e'^U), (3.33) 

we write p.ip as 

b°{eW,e'^U)dtW = -c'diW + b^{W,eU,eW,e'^U)diW+To{W,eU,eW,e^U) + eTi{W,eU,eW,e^U), (3.34) 
e 

where 

To{W, eU, eW, e^U) = foiW + d<^{W),eU, e{W + d<^>{W)), e^U) 

- &°(eiy, €^W)d^{W, eU, eW, e^U) + b\W, eU, eW)did^W) (3.35) 

and 

Ti(W, eU, eW, e^U)) = fi{W + d^{W),eU, e{W + d<^>{W)), e^U) . (3.36) 

Proposition 3.4. Suppose -1 < 5 < -1/2, eo > 0, .s e Nq, R > 0, Ki < Kq/ t > 
2Ki/Csob, R > 16t + R, k > 3 + s, a,w^ & Hi,, suppa C Br, £ e Hi,. Let ul^ , 

000 O 

dtU^J and be the initial data constructed in Provosition \3. li which, by choosing Cq < 1 small enough, 

o o 

satisfies 



(edtul^,diul^ -dlSidiA~^p,0,a,wl) < , and ||u*^|Ufc+i < 

\ O O ^00/ //fc O 5 



12 



for all e € (0, eo]- Then there exists a T > independent of t ^ (0, e^], and maps 
W, = {n%, W}% K', "e, wlf e XT„s,k,s-i < e < eo 

such that 

(i) T,>T forO<e< eo, 

(ii) We is the unique solution to (|3.34p with initial data 

\ O O Q o o J 

(Hi) 

o-l.e -'-'6-1, e 

and 

max{||eu:^(t)||i~,||ea,(i)||L-,||ew'(t)||L-} < 2Xo 
for all (i,e) G [0,T] X (0,eo], 

(iv) if 

limsup ||VFe(t)||vi/i,oo <oo, 

and 

sup { II eu:^' (i) II , II ea. (Oik- , II e^;' (0 II } < 2ifo , 

0<t<Tj 

i/ien t/ie solution VFe(t) can &e uniquely extended for some time T* > T^, 

(v) for any time which is strictly less than the maximal existence time and for which 

sup {\\ev^^{t)U^,\\eaS)\\L^AWmL^}<'^Ko 

0<t<T, 

holds, the support of satisfies 

snpp a,{t) C Bji^ yte[0,f,], 

where R^^ := IGsupg^j^j^ ||?i;f (t)||L~ + R, 

(vi) suppae(t) C Bff for all {t,e) G [0,T] x (0, eo], 

(vii) dtul^ = e-^ui^, and djul^ = Wli + Sldidi<i>{ae), where v-l^ = v-l^ + e~^Su'^ , 

' ' o 

(via) the triple {u*-' , , } determines, via the formulas (|1.7p . (II. 9p . (j2.ip . and (j2.2p . a solutions to 
the full Einstein- Euler system (|l.ip in the harmonic gauge p.4[) on the spacetime region — 
X [0,T], and 

(ix) the conclusions (vii)-(viii) continue to hold on any region of the form —S? [0, 2^] provided 
suppa,(t) C Bfj for allO<t<fe. 
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Proof, (i)-(iv): Given the initial data satisfying 
(edtui' , diul' - Sl6idiA-^p, 0, a, w 

\ O O o o 



< Ki , and ||u£^ ||^fc+i < 



CoCsoh 



for all e S (0, eo], it is not difficult using the inequalities p.23p and p.24p . and Lemmas 13.2113.31 and lA.7l 
to verify that ||PV£(0)||^fc ^ < Ki, \\dtWe{0)\\ f^k-i < 1, and the evolution equation (|3.34[) satisfies the 

conditions (|B.3P - (|B.5|) . Therefore, it follows directly from Theorem IB.ll that there exists a time T > 
independent of e e (0, eq] such that |j W^e(i)llff|_^ <2Ki < 2Ka/{^Csob), and ||5tVF,(i)||^k-i < 1 for 
all < i < r. This proves (i)-(iii). Statement (iv) also follows directly from Theorem lB.il 



(v)-(vi): Statement (v) follows from a slight modification of Lemma 7.2 in [15] while (vi) follows directly 
from (iii) and (v). 

(vii)-(ix): By (vi) we see that V,{t) = We{t) + d<^>{W,{t)) satisfies ([33]) for (i, e) £ [0,T] x (0,eo]. Then 
the same arguments used to prove (ii) and (iii) of Proposition 6.1 in [15] can be employed to prove the 
statements (vii)-(ix) of this Proposition. □ 

3.3 Zeroth order equation 

In order to discuss equations satisfied by the zeroth and higher order expansions, we will first introduce 
some notation. To begin, we define 

p 1 p p 01 p 

u-(;7, {/,...,[/), w = {w,w,...,w), 

p 1 p p 1 p 

X^iX,X,...,X), Y ^{Yi,Yi,...,Yi), 

and let 

T^{U, W) = To{W, eU, eW, e^U) + Ti{W, eU, eW, e^U) , 
B^{U, W, Y) = b\W, eU, eW, e^U)Yi , 

and 

B°{U,W,X) ^ h^\(?\J,tW)X . 

Proposition 3.5. Suppose £ > 3/2, R > 0, ~1 < rj < —1/2. Then there exists an eo > such that the 
maps 

B, : BRiH'^) X Br{H'^_,^,) x Br{H'^_,J —> H'^^,^,, 



B" : Br{H^) X Bii{H^_^ J X Bii{Hg_-^ J — > H^_i^^ 
are uniformly analytic for e G [0, eg] . 

Proof. The proof follows directly from Lemmas l3.21 13. 3( IA.7| and the fact that compositions of uniformly 
analytic functions are again analytic. □ 



14 



Next, we define 



P P-i P 1 dP 
^ ^ p\ deP 



i=0 



:F,{U{e),W{e), 



and 



where 



p p-i p p I dP 

B( U , W, Y) = -— ^^^B,{U{e),W{e), r(e)), 



P P-2 p-1 p I dP 

B( U , W,X) = -— ^^B°,{U{e),W{e),X{e)), 



p—1 P— 1 P P 

U{e) = Y,^m, VF(e) = 5^eW, X(e) = ^e«l, and Y{e) = ^e''Y. 

q=0 g=0 g=0 q=0 

Proposition 3.6. Suppose I > 3/2, R> 0, —l<r]< —1/2. Then there exists an cq > such that the 
maps 

B° : (Bn{H'^) x (if^)^'-^') x (^^(iJ^.i.J x (iJ^.^.J^'^) x (if^.^ JJ^_i,, 
are uniformly analytic for eG [0, eo]. Moreover, there exists uniformly analytic maps 



1 p—1 1 p p p 

i/iai are linear in the variables U,. . . , U , W, . . . , W, X, . . . ,X, Y,. . . ,Y, and 



P —1 

^,([/(e),VF(e))-^e4(U,W) 

9=0 



p p—1 p 
= :!^K,e(U,W), 



eP+1 



p —1 

B,{U{e), W{e), Y{e)) - £ e4( U , W, Y) 
g=0 



P p—1 P P 

= B7j,,(e, U,W,Y), 



and 



eP+i 



B0(i7(e), VF(e), X(e)) - ^'B^iV, W, X) 

9=0 



P p— 2 p—1 p 

= BVe(e, U,W,X). 
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Proof. The proof follows immediately from the Taylor expansions for J-'^, B^, and which are uniformly 
analytic by Proposition 13. 51 □ 



We note that from the definition of the above maps, it is clear that 

p ^ p ^ P P^l P p P P^2 p— 1 p— 1 

B = {W)Yi + B{V ,W,Y) and B° = X + S°( U , W , X ), (3.37) 

where 



b^{W) —b^iW, 0,0,0) and B = = . (3.38) 
With our notation fixed, we are now ready to define the zeroth order equations: 

dtW = b^{W)diW + T{W) + c^di^, (3.39) 
c^diW = , (3.40) 

W{0) = WMLo- (3.41) 

We showed in [15] that these equation are equivalent to the Poisson-Euler equations of Newtonian gravity. 
To see this, wc first note that the Poisson-Eulcr-Makino system p.l0p - (|1.12p is (non-local) symmetric 
hyperbolic, and thus we can use the results of Appendix iBl to obtain local existence of solutions. 

Proposition 3.7. Let k, s, 5, a, and w be as in Proposition \ 3.4\ Then there exists a maximal time 

o o 

Tq^ > and a unique solution 

^, e C°([0, T„''),Hti) n c\[o, To),H^:l) , 
I e C°{[Q, T^'), H^+') n C\[Q, T^'), H^+') , dtl> e C0([0, T^'),H^s-l) 

to (fmi)) - ([rT^ satisfying a(0) = a and w^(0) = . Moreover, 

o o 

Or _ £ ,^ o S. -1 A -1/0 Or- 

and 

suppa(i)ci3fl(t) vte[o,ro^O, 

where R{t) = R + t suPq<s<( ||u'''^(s)||2;,oo . 

Proof. From the weighted calculus inequalities of Appendix [X] (see also Appendix A of [15]), the Poisson- 
Euler-Makino system (|1.10p - (|1.12p satisfies the conditions required by Theorem lB.il Therefore all of the 

statements except for the estimate on the support of a{t) follow from this theorem. To prove the estimate 

on the support, we note that w' e C^{[0,Tf}''),Cl{R)) by the Sobolev inequahty ([X^ . Therefore we 





can integrate the differential equation dx^/dt — w^{t,x) to get a flow il^l{x) that is defined for all 

{t, x) e [0, To) X and satii 
equation (|1.10p implies that 



{t,x) G [0, To) X and satisfies -00 = Ir^- For each x £ M'^, define a^{t) = a{t,iljt{x)). The evolution 



^a-(t) + ^djw'{x,M=^))^''{t) = 0. 
at 2n 

By assumption, a^(0) — a{0, a:) = for all x G := \ Bji, and thus 



^^{t) = Mt,Mx)) = foi:all{t,x)e[0,T^^)xER (3.42) 
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by the above differential equation. Moreover, 

\ipt{x) - x\ < / \dstps{x)\ ds < / \w\x,tps{x))\ds <t sup |lw-'^(s)|| 



and hence it foUows from (|3.42[) that suppa(t) C Bj^^i^, where Rit) = R + t supQ<;^<f ||ti;^(s)||ioo. □ 

Using this local existence theorem, the next proposition follows by straightforward computation. 



Proposition 3.8. Let {a{t),w^ {t),^{t)} be the solution to the Makino-Euler-Poisson equations 
()1.12|) from Proposition \S.7\ and define 

W{t) = (0,-dlSil>{t),0,a{t),S}w\t)f e Xt„",.,m-i, '^nd = (^^^"(t), 0, 0, Of , 

where 



I 

Then {W{t),uj{t)} defines a unique solution to the initial value problem p.39p - p.4ip on the time interval 
0<t< T^'. 

4 First order expansion 

By Proposition 13. 11 the initial data li*^ is analytic in e and there exists a convergent expansion in H^'^^ 
for u^J of the form u*-' — X^^o s^'^*"' for < e < eo- Consequently, U can be expanded as J7 = X^^o 

q g . 

where U — (0, 0, , 0, 0)"^. Moreover, by Lemma [3T2l and the inequality (|3.24p . we can expand VFe(O) as 



W,{Q) = Y.'''W (4.1) 

q=0 



with the sum converging in Hg_^ ^ uniformly for < e < ep. 

2 

We define the second order remainder by 



W,^W + e{iJ + W,)+e^Z,, (4.2) 



1 

with the first order expansion term satisfying 



1„1 1,1 °1 °i loiO 101 101 

bJ^dtW, = -c^diW, + b^diW, + b^djuj + B(U, W, Y) - ^^(W, X) + W), 



e 

1 



(4.3) 



where 



W,{0) = W -(1j{Q), (4.4) 



b° = b°{0,eW)., 



00000 00 

V = U, W = W, X^dtW, Y = diW, 



17 



and 



W = (T^,i + V7,), X = {dtW,dt^). 



Observe that 

by Proposition [ajl Substituting (l42l) in (1X341) yields 



where 



1 



and 



(4.5) 
(4.6) 



Using we then find that satisfies 



I _ 1^,(0) -W^(0)-e(^,(0) + cL(0)) 

^t\S>) — 2 



where 



1 



9=0 

/- x-^ 9,9-1 9 9 X . ^ 9 9-1 9 

(i?. - ^ e«i?( U , W, Y)) + (.F. - ^ e«^( U , W 



and 



Letting 



9=0 



1 



9=0 



i5=(lJ,U), Y = (9/1^, a/i + 



x = (9tW^,ati + atVK,), 



it follows from Proposition 13. 61 that 

2 1 1 i 1 2 



1 2 



(4.7) 

(4.8) 
(4.9) 



(4.10) 



(4.11) 



/C, =£(U,W,X,Y,Ze) +X(e,C/,W,X,Y,eZ, 

2 

for analytic maps £ and Al with L linear in . 

As we shall see in Theorem 14.21 when the initial data is chosen such that \\dfWey^j\\ 
bounded as e \ 0, the e dependence can be removed from the first order expansion coefficient. This is 



k~2 remains 
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accomplished by replacing (|4.3p - (l4.4p with a related, but different e independent version. To describe this 
system, we let 

and define projection operators by 

Ili{W) = (l^i) and Ilj{W) = (W'j) . 
Then the system that replaces (I4.3p - (|4.4p is: 

dtW = b'diWe + b'di^ + B{IJ,W,Y) - B°{W.X)+J'{V,W) (4.12) 

W{0) -il}{0) , (4.13) 







where 



c5= (c5/^9,f^^^ 0,0,0)^, (4.14) 
Aw4 = -d''llj{B{lJ,W,Y) - SO(W,X) +i^(U, W)) , (4.15) 



and 



2 iOlO 101 101 

An = -n4(B(u, w, Y) - s°(w,x) + j^(u, w)) . (4.16) 

Existence of solutions to the initial value problem (I4.12p - (l4.13p is covered by the following Proposition. 

Proposition 4.1. Let 5, k, s, Ki, R, R, andr he as in Proposition \3.4\ Tq^ be as in Proposition \3. 7\ and 

suppose Tq < Tq^ . If s andr are chosen so that s>l, and 16t > max{32iiri, Tq supg<j<j.^ sup ||w^(i)||L°°}, 
then there exists a map 



1 



W e XTQ,s-i,k-i,& 

such that W{t) is the unique solution to the initial value problem (j4.12p - (l4.13p . and 

suppp(t)cB^ forQ<t<TQ, 

where p = p-fj77^+i)y^^^" "'^Q^- Moreover, if the initial data satisfies c^diW{0) ~ 0, then 

T ^ 2 2 

c diW{t) = for < t < To, and uj, e XTo,s~i,k-i,s-i- 
Proof. By construction, we have 

Wo - i(0) e H^,ZI ■ (4.17) 

Next, we observe that the map 

H',+' X {Hl.xHtD X Htl X Htl 3 (U,W,X, Y) 

I — *n4(i?(U,W,Y)-B°(W,X)+.F(U,W)) eH^gZl (4.18) 

is analytic for ^ > 3/2 + 1, which follows directly from the weighted estimates of Appendix [A] (see also 
Appendix A of [15]). It therefore follows that the system (|4.12p - (l4.16p satisfies all the hypotheses of 
Theorem lB.il Thus, there exists a unique solution 

W e XT„,,-i,k-i,5-i (4.19) 
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2 

satisfying the initial value problem (|4.12p - (|4.13|l . Furthermore, from p.28|) - (|4.18p . it is clear that oji = 

2 1 

diVt G XTa,s-i,k,s-i- Note that we have used the linearity of the system (|4.12p - (|4.16p in W to conclude 
that the solution can be continued as long as the coefficients are well defined, which is the case for 

< i < To < To^- 

By assumption, the initial data satisfies 

c^diW{Q) = 0, (4.20) 

while from Proposition 13.81 we have that 

W^'^ (t) = Wi'^ (t) ^ 5u'^ {t) =0, (4.2f ) 

and hence 

u^Ht)^0, ui'^t)=SlSidii>{t), and u'^t) = 0. (4.22) 



From this it follows that has a block diagonal structure of the form 
and consequently 

fe^a/i^O, U4{b'diW)=0, and nj{b^diW)=0. (4.23) 

1 

Next, a straight forward calculation using (|4?T2l) . ((4T4)) - ((4T5)) . and (|4?23l) shows that dt{c'diW) = 0, 
and hence 

c^diW(t) =0 for < i < To 
by (I4.20p . By the definition of the c^, this is equivalent to (since (5 < 0) 

W^4(i)=0 and d'Wi{t) = 0. (4.24) 
A short calculation using (I4.f 2p and (|4.24p then shows that 







dtSu'^ = Cu^'^ = 5\5idt<i> . (4.25) 
However, (5u'-' (0) = (see Proposition [SH]), and so integrating (|4.25p yields 

5u'^ =5\5i{l{t)-l{Q)), (4.26) 



and 



Also by (|4.24p . we have that 



while 



where 



u'^ [t) = u'J + 5u'^ (t) = SlSi^it) . (4.27) 

i4»J(i) = i4*J(t) + W4'Ht) = 6l6il>{t), (4.28) 

i/^(t) = i/J"(t) + SlSidj^t) + Wj'^it), (4.29) 

A$ = p. (4.30) 
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We remark that in obtaining (|4.30p . we have used supp p{t) C Bfj for < t < To, which foUows from the 
1 



(4.31) 



definition of p and Proposition 13. 71 

Using (|43T|) . gJll), gJll), g^, gJll), and ((4:29)) together, we can write (|4T^ as 

Moreover, it follows from the evolution equation ()4.12|) that 



We also note that 
by Proposition (|3.8p . and hence 



dtd[jWi]'^ = 



(4.32) 
(4.33) 
(4.34) 



by (11321). However, d[jWi]{0) = by PropositionO and thus we get from ((iTM)) that d[jWif{t) = 0. 
This combined with (|4.24p shows that (since (5 < 0) 



W^(t) = 0, 



and hence 



(4.35) 



(4.36) 



Using (USB), (11221), (1123), (jliSni), (11311), and the evolution equation (jH^ . a straightforward calculation 
then shows that the pair {a,!;;^} satisfy 



1 

Ota — w Oja — —oiw — w Oja — —uiw = . 
2n 2n 

9tW^ - w^diW^ d'^ a - d-^^ -w diW^ d'^a = . 

2n 2n 



Also, we observe that 



by (|OTl) and (|i351) . and that 



(S4 = 2<5^'ap9t(P^)+<51<5i9tCi>, 



2na^ 



and 



But, by pTTU)) . we have 



2nQ; 



2n-l 



{iKnin + 1))" 



(4ii:n(n + 1))" 



0/ 0/„ 1" 

aojw + w Oja 



w a/a H ojw 

2n 



2n(2n-l) /0/„0x02„_2: 



(4.37) 
(4.38) 

(4.39) 
(4.40) 



{4Kn{n+l)y 



dtp 



2n 



{'iKn{n+l)y 



' 1 , 2n-l i 0, 

Ota H d/w 

2n 



2n(2n — 1) . o r „ o 
(w Oja)" 



{4:Kn{n + l)y 



(4.41) 



(4.42) 



and therefore 



„1 n rO Or In 

Otp — Oi[w p + w p) =0 
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by gJT]), (|440)) . (|44T|1 . and (|442| . It then follows from (|430| that 
and hence 

2 ■ ■ 

by (mm), g3ni), and Proposition O □ 

Theorem 4.2. Lei 5, /c, s, Ki, R, R, t, T, and Wc{t) he as in Proposition \S.4[ {W{t),LLi{t)} as in 
Proposition\3M 

as in Proposition \3. 7| and suppose Tq < Tq^ . If s and r are chosen so that s >2, 
and 16t > niax{32iiri, Tq supq^^^j;^ sup ||ti;^(t)||L°= }, then for eg > small enough, 

(i) there exist constants K2, K3 such that the solution We{t) (0 < e < eg) exists on the interval [0,T^), 
where 

f. = nrin{To,i-ln(^^ 

and obeys the bounds 

sup ma.x{\\eu4t)\\L^,\\ea,{t)\\,\\ewl{t)\\L--} < 2Ko , 
o<t<f, 

sup \\We{t)\\w^.-=- < 00, snpp pe{t) C Bj^, 

0<t<T, 

(a) and there exists maps 

1 

We £ XTo,s-i,k~i,s~i < e < eo, 
such that We is the unique solution to the initial value problem (|4.3p - (|4.4p . and 

\\We{t) - W{t) ~ t{h{t) + We{t)) < \\We{t) - W (t) - e(i(i) + W^, (t)) 1 1 ^.-2 ^ < e 

for all {t,e) £ [0,f,) x (0,eo]. 
(Hi) Moreover, ifWe{0) satisfies \\d^We{0)\\^k-2 < 1 for < e < eo, then 

\\We{t) - W{t) - e(i(t) + W{t))\\Hk-2 < \\We{t) - Wit) - e(i(i) + Wit)) ||^.-2^ < e'^^'e^ 

1 

for all it,e) £ [0,Te) x (0,eo], where W £ XTo,s-i,k-i,s-i is the unique solution to the initial value 
problem 

Proof (i)-(ii): Fix < niin{r,To}, and let 

Ci= sup \\Wit)\\H^ + sup \\dtW{t)\\j^.-i, 

0<t<T, " ^ 0<i<T. 



C2 = 

and 



sup ||i(i)||Hf , + sup \\dtujit)\\fjk-i, 

0<t<T, " ^ 0<t<T, ■'-i 



C3 = ||tJ^-i(0)|| . 

U o — 1 
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Since 



1 

and We satisfies the linear equation (|4.3p . it follows from the energy estimates derived in the proof of 
Theorem IB . 1 1 that there exists a constant K2 = /^2(C'i, C2, i^o/(\/eo^Sob)) such that 

\\W,{t)\\H^-^ < e^^^*C3 + i^2 V(t,e) e [0,r,] x (0,eo]. (4.43) 



Next, we observe that 



00 d — 1 , e 

2 



< Ko + eCsob[4ZemH^:^^ + \\WS)\\„^:^^+C2] , (4.44) 
\\W,it)\\w^.^ <Csoh[e^\\km„^-- +4WS)\\„>'-^ +eC2+Ci], (4.45) 



and 



\\W,{t)-W,{t)\\w^.^ < Csob[e'||^.(i)|lH,^_-2^ +e||W,(t)||^._-i^ +eC2]. (4.46) 

22 2 
Setting Z^{t) = eZ^{t), we note that by construction there exists a constant C5 such that ||Z£(0)||^fc-2 

j-r-^ 2 

< eCi. Moreover, from the error equation (|4.8p . it is clear that Z satisfies an equation to which Theorem 

2 

IB. II apphes. Therefore, for any K3 > eC^ (0 < e < eg) there exists constants Ki,K^ such that Z[t) 
satisfies an estimate of the form 

\\Z{t)\\<e{e'^^'[Ci^^Kz]-K5)<K^ for < i < f , (4.47) 

where 

f = min(T.,-Llnf t^'f, ]]. (4.48) 

Statements (i) and (ii) now follow directly from Propositions 13.41 and 13.71 and the estimates (|3.25p . 
(H^ - fOS]) . (1133), and ((i^ . provided eq is chosen small enough. 

(iii): To prove statement (iii), we first observe that it follows from the evolution equation (|3.34p that the 

1 

condition ||9j^T/Fe(0)||^fc-2 < 1 for < e < eo is equivalent to the condition c'diW{0) — 0. Then replacing 

1 2 1 2 2 

We{t), and Ze{t) in ()4.2|1 with W{t), and u; + Z^it), respectively, it is not difficult using Proposition 14. II 

2 

to show that the new error term Z^{t) will satisfy the same type of estimate as above. We emphasize 

1 2 1 

that the key property used to make this replacement is that W{t) and u){t) satisfy c^diW{t) = and 

2 1 1 

uj & XTo,s-i,k-i,s-i- The proof of statement (iii) now follows as we are able to replace W^(t) with W{t) 

everywhere in the above estimates. □ 

5 Higher order expansions and convergence 

1 

Theorem 5.1. Let S, k, s, Ki, R, R, and W:(t) be as in Proposition \3.4\ {W{t),uj{t)} as in Proposition 

1 

as in Proposition \3. 7[ W ^ [t) and r as in Theorem \4-2\ and suppose To < Tf}'. If s > 3, then 
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for Co small enough, there exists an infinite sequence of maps 

q 

We e XTo,s-2,k-2,s-i q e Z>2 

such that 

(i) each W{t) satisfies a linear (non-local) symmetric hyperbolic system with initial data Vl^e(O) = W 

]^ r r 

and coefficients depending on e, W , oj, U for < r < q, and W ^ for 1 < r < q — 1, 

\\W ,{t)\\HU-^ ^ e\\dtW ,{t)\\HU-. < \\W,{t)\\j^.-2^+e\\dtW,\\j^.-3^ < 1 
for all {t,e,q) G [0, Tq) x (0, Co] x ^>2, and 

(Hi) 

^ oo 

W,{t)^W{t) + e{i!u{t) + W,)+J2^'^^'^*^ (t,e) G [0,ro) X (0,eo], (5.1) 

9=0 

where the sum converges uniformly in {[0,Tq)] HgZi e) and {[0,Tq); H'^^'^) . 
(iv) Moreover, if s — 2> p>l, and the initial data is chosen so that 

l|5r'W^.(0)||^'.-(,+i) <1 q=l,2,...,p, 

then there exists e-independent maps 

« 9+1 

W & XTo,s-q.k-q,5-i and lo E XT„,s^q,k-q,S-i q=l,2,...,p 

such that 

(iv.a) each W satisfies a e-independent linear (non-local) symmetric hyperbolic system with coeffi- 

r J, r 

dents depending only on U for < r < q, co for < r < q + 1, and W for < r < q — I, 
and 

1 q 1 

(iv.b) the terms W ^ in the sum (IS.ip can be replaced by lu + W for 1 < q < p with the sum converging 
uniformly C°{[0,To), H^Z[l)^^'') and C°{[0,To); H''-^'i+^'^) . 

Proof. The proof of this Theorem follows from a straightforward adaptation of the proof of Theorem 3 
in [21]. We will only sketch the details. 

Following Schochet [21] (see also [11]), we consider the following iteration: 

rn m+1 1 m+1 m m+1 m+1 m 

b°{Z,)dt = -c'di Z, + b\Z,)di Z, +C{Z,) + eM{Z,), (5.2) 
Ze (0) = e^-'W, (5.3) 

9=2 

where 

mo -1 1 m m m rn 

Zi=0, We^W -^e{uj + We) + e^Ze, b^iZe) = b^{We,eU,eWe,e^U), 

771 IP' m+1 2lilm+l m lilm 

b°{Z,) = b°{e^U,eW,), £( Ze ) =£(U,W,X,Y, ), and M{Z,) ^ M(e,U,W ,X,Y , Z,) . 
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Using the energy estimates of Theorem IB.ll and the weighted Sobolev estimates in Appendix [X] (see 
also [15]), it is clear the arguments of Schochet can be generalized to show that 

m m 

\\Z{t)\\^.-2^ + e\\dtZ{t)\\^.-3^ < 1, (5.4) 

and 

ra-\-2 m+1 m+1 m m+2 

II Z (t)- Z (t)ll^. 3 <e|| Z (<)-Z(t)||^. 3 +e"|| II . 3 (5.5) 

— 1, e — 1, e (J — 1,6 

for all (i, e) G [0,To) x (0,eo]. Therefore by (|4.ip . (|5.4p . (|5.5p . and the uniqueness of solutions to the 

1 1 m 

evolution equation (|5.2|) . we see that for ep small enough the sequence W{t) + e(We{t) + u!{t)^ + e'^Z^{t) 
converges in C^i[0,To),H^I^J to We{t) for each e e (0,eo]. Therefore, defining 

m+1 m 

w\it) = 

we have that 

1 1 ^ q 

W,{t) = W{t) + e{W,{t) + uj{t)) + e''We{t) 

q=2 

with the sum converging in C^dO, To), ^^^Zi J for each (e e (0,eo]- Moreover, because of the inequality 
(j3.25p . it follows that the sum converges uniformly in C"([0, Tq), i?''-^) for e e (0,eo]. This completes 
the proof of statements (i)-(iii). The proof of statement (iv) also follows easily from the arguments used 
in the proof of Theorem 3 in [21]. □ 

Remark 5.2. The equations satisfied by the W from part (iv) of Theorem ()5.ip are: 

9 9 „ i 0-1 9 9-1 

dtW ^b\W)diW + b\W)diCo + B{\J Y ) - dtuj 

in/9-2 9-1 9-1, 9 ,9-1 9 , , n+i 

-s°(u, w, X) +j='(u,w) +c^a/^i , 

c^djW = , 
W{Q) = W -uj{Q) , 

where 

9 q q Oil n 1 
U= ([/,...,[/, W ^{W,uj + W,...,uj + W), 

^^{dtW,dtlu + dtW,...,dtU} + dtW), and Y = {diW ,dilu + diW , . . . ,diuj + diW). 

6 The first post-Newtonian expansion 

We are now ready to prove the main theorem that guarantees the existence of a large class of solutions 
to the Einstcin-Euler equations that can be expanded to the first post-Newtonian order. 

Proof of Theorem Using the harmonic equations 

ea,u44 ^ „aju4/^ ^Q^-IA ^ _Qj^iJ^ (6.1) 

we can write the constraint equations p. lip as 

Au^*^- = 5lp - S'^jdLedtu" + e\cif{eu'\didju'^ ,edidtu'^'') 



25 



where Qq^ (2/1, 2/2, 2/3 ) is bilinear in yi and (2/2, ys), '5i''(yi, 2/2, J/s) is quadratic in 2/2,2/3, and the maps Qf" 
[v = 0, 1, 2) are analytic in all their variables for e^u^^ G V. We can also write the fiTL-components of the 
reduced Einstein equations p.l3p as 



,2-KL 



£2(1 



(6.2) 



where Q^^(2/i, 2/2, 2/3) is quadratic in (2/2,2/3), 

<9f^ = Qf^(eu'^e^u'^w,ew)a2 + Qfi(euy\e2-y,^,ew)i«^w^^ 

and all of the maps {y — 0, 1, 2, 3) are analytic in their arguments for (^u G V. 

We now take 

as the prescribed initial data, and solve the non-linear elliptic system 
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Au^^ = A^^ := ~2e^u^^didtu^^ - e'^u^'^ Diju^^ - e^Q^^ {e^u'\ Omu'^ edtii^'^) 



(6.3) 



(6.4) 



to determine the initial data {u^-' \t=o, dtU^^t=Q} on So= {{x^ ,0) \ (x^) £ M'^}. Note that is determined 
by the fluid velocity normalization (|3.13p . which can be written as 



(6.5) 



where 7(2/1,2/2) is analytic in a neighborhood of (0,0) and /(y) = 0(|y|2) as y ^ 0. 

Using the weighted multiplication inequality (see [15], Lemma A. 8) and Lemma [A. 7| it is straightfor- 
ward to verify that there exists an eo > such that A'^ (see (|6.3p - ()6.4p ) defines an analytic map 



{e,iy,a,w',f",u'^) e (-60, eo) x Hi, x H^l x Ht, x H^^ x A^^ e H^^ , 



where 



Writing ^^Elj-^M as 



A^' = S^p + 0(e) and A^^ ^ 0(6^) as e \ 0. 



(6.6) 



it follows from (|6.6p and the invertibility of the Laplacian A : Hg H^Z^ that we can use the analytic 
version of the implicit function theorem [8] to conclude that there exists an open neighborhood U of any 
point in Hg_, x Hg_, x Hg_, x H^Z^^ and analytic maps 

{e,ii' ,a,w' ,\") e (-eo,eo) x t/ ^ u''^ G HI 



that solve equations (I6.3p - ()6.4p . Moreover, it follows from ()6.6p that 

l|uf^(0)||i^,^ V6e[0,6o], 



(6.7) 
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and hence 

\\dM'^m\Hl_,<^^ Vee[0,eo]. (6.8) 
Also, we note that by construction 

l|ftaf^(0)||^.-i <e2 Vee[0,eo]. (6.9) 
Differentiating the harmonic conditions (|6.ip with respect to t, and using (|6.7p - (|6.9p . yields 

l|5fflr(0)||^.-. < 1 p = 0,...,4, (6.10) 

and 

||9fuf (0)11^.-. <e p = 0,...,3 (6.11) 

for aU e £ [0, eq]- 

Using (|6.ip . the Euler equations (|2.26p can be written as 

Sjw = [a4(e2uy , ew)] (a^(w, e^u^ , ew)a/w + 6o(a/u'J , eajU^^)^ 

5i(w,e2u'^ew,9/u*^eatu^^,ea/u'^e2(9tu^'^)) , (6.12) 

where the maps a^, , bo, bi are analytic in all their arguments for e^u S V, and a"'^(0, 0) = I , 

a^(0, 0,0,0) = 0, bo{yi,y2) is linear, and 6^(yi, y2, ?/3, 2/4, ^5, 2/6, y?) is h ncar i n (2/4,2/5,2/6,2/7) and sat- 
isfies 6^(0,0,0,2/4,2/5,2/6,2/7) = 0. Then differentiating (|6.ip . (|6.2p . and (|6.12p with respect to t while 
using (|6.7p - (|6.1ip shows that 

|iafuf^(0)||^._-p <1 p = 3,4, (6.13) 

||9fae(0)||^._-p < 1 p = 0,...,3, (6.14) 

||«(0)||^.-p < 1 p = 0,...,3, (6.15) 

and 

\\edy\\„.-.<l (6.16) 

for aU e e [0, eo]. We then find from the definition of W^, the estimates (IST7)) - (|S11|) . and (|0^ - (p7TB)) . 
that 

ll'9t^^e(0)||^^-P < 1 for j3 = 0,1,2,3 and < e < eo. (6.17) 

Next, we observe that 

by (p:^ and (0271), while for any < ^ < fc, 

mm < \\vM\Hi_,^^ = \\We{t) + d^{w,m\Hi_,^, < \\wmHi_,_^ (6.19) 

by (1221), (12211), and (P3T|) . The proof of Theorem O now follows directly from Theorem O and the 
estimates ([CT7|) - ([Sjg| . □ 
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7 Discussion 



In this article, we have estabhshed the existence of a large class of dynamical solutions to the Einstein- 
Euler equations that have a first post-Newtonian expansion. Although this is an improvement over 
existing rigorous results [15,19], which only cover the Newtonian limit situation (i.e. the "zeroth" post- 
Newtonian expansion), the results of this paper are almost certainly not optimal. In general, one expects 
that with a suitable gauge choice, it should be possible to generate post-Newtonian expansions to at least 
the 2.5 post-Newtonian order after which there are indications that the post-Newtonian expansions will 
break down. For a lucid discussion of this phenomenon see [17]. 

As remarked in [17], the choice of harmonic gauge may be the reason for not being able to reach 
the 2.5 post-Newtonian order. At the formal level, there exist other gauges that perform better than 
the harmonic gauge for the post-Newtonian expansions. However, it remains to be seen if these other 
gauges are compatible with the singular hyperbolic energy estimates that are guaranteed to arise in the 
dynamical setting. We are presently investigating this problem. 

From the proof of Theorem 1 1.1 1 and the paper [15], it is clear that conditions of the form 

\\d?W,{0)\\H^-p < 1 as e \ (7.1) 

on the initial data play a crucial role in generating the post-Newtonian expansions. This leads to the 
question of what happens when one considers initial data that does not satisfy (|7.ip for any p G Z > 0. 
In [16], we address this question for the situation where 

hmsup ||9tM^e(0)||^fe-i = oo. 

There we find that a Newtonian description is still appropriate for the motion of the matter, but the 
gravitational field no longer vanishes in the limit e \ 0. Instead, there exists high frequency gravitational 
radiation that is not small at the e° order, and this will necessarily affect the higher order expansions. 
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A Weighted calculus inequalities 

In this section, we prove additional weighted calculus inequalities that are similar in spirit to those in 
Appendix A of [15]. We first recall from [15] the definition of the weighted Sobolev spaces. Let y be a 
finite dimensional vector space with inner product (•]•) and corresponding norm | • ]. For u G L'^^^{W^ , V), 
1 < p < oo, (5 G K, and e G K>o, the weighted norm of u is defined by 



^^s "/p^ii^^ if 1 < p < oo 



(A.l) 



if p = oo 



where a^ix) 



1, 



1 + — jea;]^. The weighted Sobolev norms are then defined by 



(E ii^"<^,J'^' ifi<p< 

if p = oo 



\a\<k 



(A.2) 



I |a|<fc 
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where k E Nq, a — (ai, . . . , q;„) E Nq is a multi-index and D" — d"^ . . . d"" . Here 

where (x^, . . . , a:") are the standard Cartesian coordinates on M". The weighted Sobolev spaces are then 
defined as 

W^f ^{ue Wt\R\ V) I \\u\\^... < ex. } . 

We note that W^'^ are the standard Sobolev spaces, and for e > the Wgf are equivalent to the radially 
weighted Sobolev spaces [1,7]. For p = 2, we use the alternate notation Hg^^ :— Wg'^. The spaces 
and Hg ^ are Hilbert spaces with inner products 

{u\v)l2 / {u\v)a-^'--d'^x , (A.3) 

and 

|a|<fc 

respectively. When e = 1, we will also use the notation Wg''' = Wgf and = i?^^. 
Lemma A.l. Suppose eg > 0, (5i > max{(52 + (^s, 5/^ + (Js}, f/ien 

for all € e [0, eo],ue Lf;^^ n ilf^^,, and v G i,?^^, H i/f^^,. 
Proof. This follows directly from the inequality 

\\uv\\Hk < ||li||Loo||w||^fc + \\Du\\Hk-i\\v\\L'>- 

and Lemma A. 4 of [15]. □ 

Lemma A.2. Suppose eo > 0, 6 < 0, -n/2 < A < -n/2 + I, X > 5, k > n/2, and f e C^{^^ x 
M^, M^^^"'^"') with /(O, 0) = 0. Then there exists a polynomial p{yi,y2,y3) such that 

\\f{u,w)v\\jfk^^ < \\f\\c^p{\\u\\„k,\\w\\„k^^,\\v\\„kJ\\v\\Hi^^ 

for all e £ [0, eq], u & and w,v £ Hg^. 

Proof. Since S < X, it follows from Lemma [A. II that 

\\f{u,w)v\\H.^^<\\f{u,w)U^\\v\\H^^^ + {\\D{f{u,w))\\j^.^^^^^^ 

Using Lemma A. 9 of [15], we can write the above inequality as 



+4AhI^ + l|u'||ffj_-i)||f|lL- . (A.5) 
But k > n/2 and X< 5 <{) implies that 

hlU- < , IIw'IIl- < , M\l-- <\\v\\LT,'^M\m , (A.6) 
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and 



< Mm (A.7) 



by equation A. 24 and Lemma A.7 of [15], while 

\\Du\\jj.-i^+e\\uh2^<\\u\\H^^ (A.8) 

follows from Lemma A. 11 of [15] since —n/2 < A < — n/2 + 1. The proof now follows directly from the 
inequalities (|A.5|) - (|A.8p . □ 

Lemma A.3. Suppose > 0, (5 < 0, -n/2 < A < -n/2 + I, X > S, k > n/2 + 1, and f e C^'(K^ x 
]^Af^ jyjMxM-j ^j^/j /(0,0) — 0. Then there exists a polynomial p{y 1,1/2) such that 

\\[D°'J{u,w)]v\\l2_^^^ < \\f\\c^pi\\u\\H^^,\\w\\Hl){\\u\\Hk + \\w\\HkJ\\v\\^k-i^ 

for all e S [0, eo] , 1 < |a| < k, u e H^, w G Hg^, and v G HgZi ^■ 

Proof. The proof follows directly from Lemma A. 9 of [15] and the inequalities ()A.5P - (jA.8p . □ 

Lemma A. 4. Suppose eo > 0, S < 0, —n/2 < A < —n/2 + 1, X > S, and k > n/2. Then there exists a 
constant C > such that 

||"lW2||fl-!; < C||Mi||^j||M2|liJj , 
\\uiVi\\hI^ < C||Mi||^fc||wi||^fc^ , 

and 

\\viV2\\hI^ < C\\vi\\hiJ\v2\\hI^ 
for all Wi, 1*2 G H^, Vi,V2 G Hg^, and e G [0, Eq] . 

Proof. The proof follows immediately from Lemma [A. II and the inequalities (jA.6|) - (jA.8p . □ 

We now recall the definition of analytic maps between Banach spaces. 

Definition A. 5. Suppose Y and Z are Banach spaces, U C y is an open set, and Cj{Y, Z) is the space 
of continuous, j-multilinear maps from Y to Z with norm 

\\F\\cj(Y,z) = sup {||F(ui,M2, ...,Uj)\\z\ Uj G U and sup{||ui||y, ||u2||y, . . . , Hualk} < l}- 

Then a map / : U — > Z is analytic in U, if for each uq ^ U there exists a p > 0, and a sequence of maps 
multilinear maps fj G CjiY, Z) such that 

X! \\l3\\c,(Y.Z)P^ < 00, 

and 

00 

f(u)^^fj{u-uo,...,u-uo) (A. 9) 

j=o 

for all u G U satisfying \\u — uo\\y < p. The set of all analytic functions in U will be denoted C^lU, Z). 

In addition to analytic maps, we will need analytic maps that are uniformly analytic on the H^^ 
spaces as e varies. 
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Definition A. 6. Suppose R > 0, Y, Z are Banach spaces, and V dY is open. Then a sequence a maps 
fe : Br{H^^) X V ^ e ^ ^ ^^^^ called uniformly analytic for e £ [0, eo], if 

(i) /, e C^iBniH^l^ X V; H^l^ x Z) for < e < eo, and 

(ii) for each vq € V there exists constants p, cj > 0, and a sequence of maps multilinear maps /| € 
^li^k.e ^ ^'^sle X ^) such that 

00 

Y,cjip+Ry < 00, 

j=0 

and 

00 

«) = X! ''^ - I'o, ■ ■ ■ , M, I' - Wo) < e < eo 

i=o 

for all {u,v) G Hg^ ^ x V satisfying ||u||^fci < R, and \\v — vq\\y < p- 

The next lemma shows how to construct a particular class of uniformly analytic functions. 

Lemma A.7. Suppose eg > 0, (5 < 0, -n/2 < A < -n/2 + l,k> n/2, F e C"^(Bhi(M) x B;v^,(M),R), 
0) = 0, and C is the e independent constant from Lemma \A.4\ Then for < e < eo, 



00 00 _j 



p=0 q=l ^ ^ 

defines a function of class ^^[B^^ i^x) x -^^^(^l e)' e) 'where Ri ~ Ri/C and R2 = R2/C 

Proof. Using Lemma IA.41 the proof follows from a slight modification of the proof of Proposition 3.6 
from [10]. □ 

We note that the above Lemma can be easily generalized to maps / € C"(i3_R(R^) xi3ij(R^^), Mmxm)- 

B Symmetric hyperbolic equations 

The hyperbolic equations that we will consider are of the form 

b°{eue,ewe,eve)dtv^ = ^c^ d^v^ + V {e,u^,w^,v^)djV^ + 'yF{€,Ue,w^,v^), (B.l) 

Ve\t=0=V„ (B.2) 
o 

where 

(i) the maps = ^^(x) and = Wi{t,x) are and valued, respectively, while the map = 
We (t, x) is R*^-valued, 

(ii) F is a (possibly non-local) map satisfying 

\\F{e,U,Wi,Vi) - F{e,U,W2,V2)\\Hl ^p,eo,kJ \\wi - ^211^^ + - V2\\h'I (B.3) 

for aU e G [0, eo], u e Bp{Hl), wi,W2,vi,V2 € Bp{HlJ, and 

\\F{e,u,w,v)\\Hi^ <pi\\u\\„e^,\\w\\„k^^,\\v\\jfk ){\\w\\„k^^ + MhiJ (B.4) 
for aU e e [0, e], u G Hi, and w,v e Hg 
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(iii) 6°, V e Cfc^(R^ X X R''',Mmxm) (j - 1, . . . , n), 

(iv) and V are symmetric, 

(v) the o' are constant symmetric matrices, and 

(vi) there exists a constant uj > such that 

>^Imxm foraU(a,6,6)eM^xE*^ xK*^ (B.5) 

Let [n/2] denote the largest integer with [n/2] < n/2, ko = [n/2] + 2, and 

s+l 

Theorem B.l. Suppose eq > 0, T > 0, s e Nq, k = ko + s, S < 0, -n/2 < A < -n/2 + 1, 

e Hi, e H^, We G ^T,s,fe,6, < e < eo, 

o 

and 

WVcWh^ <Ci, \\w^{t)\\Hk +\\dtW^{t)\\^k-l <C2, llUelln*. < C3, 

for constants Ci, C'2, C3, independent of{t,e) G [0, T) x (0, eo] . T/ien </iere exists a polynomial p{yi,y2, 1/3) 
and maps 

Ve e Xt,,s,/c,5 < e < eo, 

such that for all e G (0, eo] 

(i) Vc(t,x) is the unique solution in L°°{{0,T^), H^) n Lip((0, T^), to the initial value problem 

(EH) -(1121;, 

(a) i/hmsupj^y^ ||we||i4'i,oo < 00, then the solution can be extended (uniquely) for time T* G [T^,T), 
(iii) for any constant Ki > Ci, 

1C2 



\W{t)\\HU^^<e^^{K2{l+i)t) 
for all {t,e) G [0,T) x (0, eo], where 



K2(l+"f) 



K2{l+7) 



K2 :=p(C3,C2,i^i), 

and 

. 1 ^ 7^1^2(1 +7) +7^2 

- ^= T' ^?^(TT^ I c,/^2(i + 7) + 7C2 

(tv) e\\vS)\\Hl^ < 1 V all it,e) G [0,f) x (0,eo], 

(iv) and if \\c^ djV,\\ j^k-i < e, then \\dtv^{t)\\jjk^i < 1 for all {t,e) G [0,1") x (0,eo]. 
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Proof. We will only prove statements (iii)-(v) as (i)-(ii) follow from a slight modification of arguments in 
Appendix B of [15jl. 

Let = D°'ve, hi = b'^{eue,ewe,eve), h{ = V {e,Ue,We,Ve), and = F{e,Ue,We,Ve). Then from the 
evolution equation (jB.ip . we find that 



Differentiating this yields 



where 



1 



= bl[D-^ {b1)-\e-'c>+b^)]d,v':+^blD'^m-'F:) . 
Energy estimates (see Lemma 7.1 in [15]) then show that 

^llbfllloA. ^ (lldivfo.lk^. + \\c + ebU^)\K\\Ll^ + 1^. J|<|L. , 
where divb^ — dtb^ + djb{, c— (c^, . . . , c"), b — (6^, . . . , 6"), and 

\a\<k 



(B.6) 

(B.7) 
(B.8) 

(B.9) 
(B.IO) 



Since b^ — 6"(eue, eWe, eWe), it follows from Lemma [A3] that 

||[i?",(fe?)-'(e-V+6i))]aXllLt (B.ll) 

'^P{<^\We\\Hi^^\\We\\Hl^^<^\\'Ve\\Hl){\W\\Hi + W^eWul^ + 1 1 ^'e || //^ J || 1 1 ff^ ^ (B.12) 

for some polynomial p{yi, 2/2, Hi)- Using this estimate along with (|B.4[) and Lemma lA. 21 we find that 

11/^1^2 '^P{\W\\HiA\We\\Hl^A\Ve\\Hl){'^W^\\Hk^ (B.13) 

for some polynomial p(yi, ?/2, 2/2)- Combining the two estimates (|B.9I) and (|B.13[) . and summing over a 
(0 < \a\ < k) yields 



dt 



But 



\\Ve\\k,5,e <P{\W\\HiA\We\\Hl^AMHl){l\\We\\Hl^ + (1 + 7) II J 



kellffj < ||be|||fc,5,e < 1 1 1 1 CP 1 1 ^'e 1 1 //f 



(B.14) 



by (jB.sp . and so it follows from (jB.14p and Gronwall's inequality that for any constant i^i > Ci, if we 
let if2 =p(C3,C2,ifi), then 



\W{t)\\m <exp(i^2(l + 7)i) 



1C2 



1^2(1 + 7) 



1C2 



K2{l+l) 



(B.15) 



for all t such that ||fe(t)||jj| < Ki. But HwjUvi^i.oo < ||t;e||jLfA: by Lemma A. 7 of [15], and hence, by the 
continuation principle (ii), we see that 



Tf > T :^ min < T. 



1 



' ^2(1 + 7) V CiK2{l + 7) + 1C2 



In 



(B.16) 



^The only real difference is the proof of the convergence of the Galerkin approximations. For the non-local problem, 
one can use the global compact imbedding Hg C (k > £,S < rj) to obtain convergence instead of the local compact 
H^(Br) C H{BiiY (fc > ^) imbedding used in [15]. 
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Next, differentiating (|B.f p with respect to t, it is clear that dtv^ satisfies a hnear equation of the same 
structure as (|B.fp . and therefore the same estimates used to derive (jB.lSp also show that there exists 
constants K2, such that 

\\dMt)\\H^ < e^^'||atz;e(0)||^.-i + A-3 V {t,e) e [0, T,] x (0,eo]. (B.17) 
The proof now follows from the estimates (|B.15P - (|B.17p . □ 
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